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High thermal conductivity is required for high performance electronics and one advantage of car- 
bon nanomaterials as candidates for future electronics is that they have peculiar thermal conduc- 
tivites. In this letter, we address the effects of Anderson type disorder to lattice thermal conductivity 
through long carbon nanotubes (CNTs) using atomistic Green functions. The dependence of phonon 
transmission is analyzed as a function of the length of the disordered CNT, thermal conductance 
as a function of temperature is calculated for different lengths. Analysis of transmission amplitudes 
as a function of length of the disordered device shows that phonons with different energies display 
different transport regimes, i.e. quasi-ballistic, diffusive and localization regimes coexist. In the 
light of the results we discuss heating of the semiconductor CNT device in electronic applications. 



INTRODUCTION 



Thermal transport in low-dimensional systems is of 
central importance for applications. By properly control- 
ling thermal properties it is possible to enhance the de- 
vice performance. In electronic applications, high values 
of both electronic and phononic conductances are desired. 
On the other hand a low phonon conductance is required 
in order to have efficient thermoelectric energy conver- 
sion. In semiconductor materials heat transport is gov- 
erned mostly by phonons and it is previously shown that 
the thermal conduction is strongly influenced by device 
dimensions [1]. For example, it was observed recently 
that phonon conduction can be reduced without effect- 
ing the electronic transport importantly in Si nanowires 
which leads to an enhancement in the thermoelectric fig- 
ure of merit by two orders of magnitude compared to its 
bulk value [2I, H, . This decrease in thermal conduction 
is mainly due to the scattering of phonons at disordered 
surface whereas electron transfer is maintained by bulk- 
like states. Thus, surface to volume ratio is a parame- 
ter to control the transport properties of these devices. 
Graphitic allotropes have the exceptional property that 
they are one atom tick. There is a growing interest in the 
field of phononic energy transport through carbon based 
materials. It is shown that thermal conductance of nan- 
otubes can be tuned by sliding the inner shell inside the 
outer shell [5]. Disorder induced localization is believed 
to be a possible explanation of the exponential depen- 
dence of thermal resistance on the telescoping distance. 
Conductance is independent of the device length if the 
device is pristine. Real systems, on the other hand, al- 
ways include disorder which we model with a distribution 
in force constants. Recently the effect of isotopic disor- 
der on thermal conduction through nanotubes is studied 
theoretically [6j] and good agreement with experiment [t^ 
is achieved. 

In this letter, we analyze the effects of disorder on ther- 
mal transport through semiconductor CNTs. First we 
summarize the atomistic Green function method in cal- 



culating the transport properties. Then we investigate 
the dependence on length of the device and on the oper- 
ating temperature. The energy dependence of transport 
regimes are shown to coexist and their role on device per- 
formance depending on system parameters is discussed. 



METHOD 

As it is common in transport calculations, we apply 
the partitioning scheme and divide the system as left, 
center and right regions (L, C, and R respectively). In 
the harmonic approximation the vibrational Hamiltonian 
can be written as 



RC^^c^ 



(1) 



where the first term i^^. = 1 /2 (i/^ | i/" ) + 1 /2 (ix" | i^"" | ii" ) 
is the Hamiltonian of the decoupled subsystem a, the sec- 
ond and third terms stand for the coupling between the 
central region and the reservoirs. Here \u^) is the vector, 
and (u^\ is its transpose, of mass renormalized coordi- 
nates = yJraiX^ , being the i^^ degree of freedom of 
subsystem a with the mass associated to this degree 
of freedom. K^^ is the matrix element representing the 
coupling between mass renormalized coordinate i of sub- 



system Oi with j subsystem and K^^ = kfj / -^rriimj 

with k^j^ being the spring constant in direct coordinates. 
In our calculations, we make use of force constants with 
fourth nearest neighbor parameters which are shown to 
have good agreement with density functional theory cal- 
culations [8| [see Fig. [1]. 

We calculate the phonon transmission using atomistic 
Green functions. Retarded Green function for a subsys- 
tem in the absence of coupling to the reservoirs is de- 
fined with g'^^^^^^luj] = [(cJTiO+)^/-i^""]"^ with / 
being the identity matrix. Self energy due to coupling 
to reservoir a is = K^o^gria),aj^aC 

tarded function for the central region in the presence of 
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FIG. 1: Phonon dispersions of graphene (a) and CNT(10,0) 
(b). Lattice structures are shown in the upper part; unit cells 
are indicated with dashed lines. 

reservoirs reads 

= [{co + i0+)2/ - K^^ - S^'^ - S^'^] . (2) 

In our calculations, we calculate the free Green function 
of the semi- infinite reservoirs using the surface Green 
function matching method [9] and the broadenings with 
the relation = i (S^'" — E^'"). Thermal conductance 
is defined as 

CO 

^{T) = / ^ ^(^) (3) 



where T{uj) = Tr [G'^T^G'T^] is the transmission func- 
tion, /b is the Bose distribution with co being the angular 
frequency and T being the temperature. We introduce 
Anderson type disorder by a random distribution of rrii 
at the C region which transforms to a distribution of dis- 
ordered coupling matrix elements . Green functions 
of the disordered region are obtained using decimation 
techniques and the transmission functions are obtained 
by averaging over 250 sample configurations. 

RESULTS AND DISCUSSIONS: 

First of all we note that, unlike its electronic counter- 
part, phononic heat transmission takes place not only in 
a small energy window but phonons of all energy values 
contribute to the conductance. In this sense, low energy 
phonons play a special role in phonon thermal conduc- 
tion in both pristine and disordered systems as it will be 
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FIG. 2: (a) Schematics of the system under consideration. 
(10,0) CNT is placed between reservoirs which have a tem- 
perature difference dT. We include Anderson type disorder 
at the central region whose length L is varied, (b) Ther- 
mal conductance as a function of temperature T for (10,0) 
CNTs with diffrent lengths, (c) Transmission versus energy 
for CNTs with diffrent lengths L (solid curves). Black curve 
is the transmission function of pristine (10,0) CNT. Dashed 
curves are the prefactors p(a;,T) for T =25, 100, 300 K (see 
text). (Same color code is used in subplots (b) and (c)). 



discussed below. We consider (10,0) CNT which is placed 
between reservoirs made of the same material [Fig.[2fa)]. 
The reservoirs are kept at different temperatures T and 
T -\- dT. We introduce Anderson type disorder in the 
central region and analyze the dependence of the trans- 
mission function and thermal conductance to disorder for 
CNTs with different lengths. The transmission function 
for the pristine CNT reflects the fact that each phonon 
branch is contributing to transport with unit probabil- 
ity. Namely, the transmission function at energy co is 
the number of phonon branches N{ij) crossing at this 
energy each acting as an independent transmission chan- 
nel [Fig. mb)]. Increasing the length of the C-region, 
the transmission is reduced but the low energy transmis- 
sion is almost unaffected by disorder [Fig. [Sfc)]. On the 
other hand high energy phonons are extremely sensitive 
to disorder in the sense that the transmission is blocked 
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tion at different lengths depending on their energies as it 
is shown in Fig. [3^b). Different transport regimes coex- 
ist for phonons with different energies. Altough high en- 
ergy phonons experience locahzation, a pure locahzation 
regime or an exponential dependence of thermal conduc- 
tance on CNT length is not observed. 

Therefore heating of the semiconductor CNT device 
in electronic applications is dependent on the length of 
the device and the frequency distribution of generated 
phonons. Low energy phonons will be discharged effec- 
tively while localization of energetic phonons may cause 
overheating of long devices. Localization induced heating 
becomes more dramatic at high operating temperatures. 



FIG. 3: Different regimes of phonon transport through (10,0) 
CNTs. Normalized transmission values are plotted as a func- 
tion of length for various energies. At low energies (25 cm~^, 
50 cm""*^) T((jj)/N((jj) can be fit to a line (a). Intermediate 
energies like 100 cm~^ and 200 cm~^ deviate from the linear 
behavior and show 1/L behavior (a). High energy phonons 
are best described with an exponential law (b). 



to a great extend even for the shortest regions with dis- 
order. We also note that disorder induced scattering is 
more important near van Hove singularities in the density 
of states. We observe that thermal conductance drops 
as the length of the central region is increased but this 
drop is significant only at high temperatures [Fig. [2](b)]. 
As the temperature increases the difference in conduc- 
tance for different lengths increase and then saturates. 
The low temperature thermal conductance is insensitive 
to disorder because of two reasons. First, because the 
transmission of low energy phonons itself is insensitive 
to disorder. Second, the termp(cc;,T) = uo dfE^oo^T) / dT 
in the integrand of Eqn.[3]is filtering out the contribution 
of high energy phonons to low temperature conductance. 
At high temperatures p{oo) flattens and the filter effect 
diminishes. 

In order to identify the different transport regimes, we 
analyze T(cc;, L)/N{uj) for a number of u values, and as a 
function of L. Transmission amplitudes are normalized in 
this way using their pristine values N{u;) in order to en- 
able comparison. In Fig|3ja), it is shown that low energy 
transmission {cj =25, 50 cm~^) decreases linearly with 
L reflecting the fact that low energy phonons display a 
quasi-ballistic behavior. For higher energies {uj =100, 200 
cm~^) the transmission function deviates from the linear 
fit to short device transmission (see Fig. [31(a)) and a dif- 
fusive behavior is observed. On the other hand, energetic 
phonons (cj =800, 1200, 1500 cm~^) experience localiza- 



CONCLUSIONS 

In summary, we analyzed the effects of disorder on 
phonon thermal transport of semiconductor CNTs. We 
show that different transport regimes for phonons of 
different energies coexist and their relative weight on 
thermal conductance depends on temperature. Long 
wavelength phonons are transmitted effectively but short 
wavelength phonons are likely cause overheating and the 
effect localization will be more pronounced at higher op- 
erating temperatures. 
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